PHEW 


ARID LAND GEOGRAPHY 


doi;10. 12118/j. issn. 1000 - 6060. 2020.01.11 


STE BE HDX AOD 时 空 变化 及 影响 因子 的 地 理 探测 


z 1,2 
= M”, 


JN 


孙 艳 玲 ， 高 


Ke’, 


(1 天 津 师范 大 学 地 理 与 环境 科学 学 院 , 天 津 300387 ; 


陈 #)', @ 隆 ， 马 4' 
2 西安 地 球 环境 创新 研究 院 , 陕 西 ”西安 710061) 


摘 要 : 基于 MODIS 3 km AOD 遥感 数据 ,利用 空间 自 相 关 模 型 及 地 理 探测 器 对 2010 一 2016 年 京 津 
KAX AOD 的 时 空 变化 特征 及 其 影响 因子 进行 探测 。 结 果 表 明 :(1) 2010—2016 FRP EWE E 
平均 AOD 值 为 0. 83 ,其 中 天 津 市 年 均 4OD 值 为 研究 区 最 高 ,河北 省 次 之 ,北京 市 最 低 。 研 究 区 及 
各 分 区 40D 的 7a 变化 趋势 大 体 一 致 ,呈现 先 下 降 后 上 升 再 小 幅 波动 的 状态 。(2) 空间 自 相关 分 


析 表 明 京 津 费 地 区 AOD 空间 分 布 呈 


现 显著 正 相关 。 局 部 高 高 聚集 区 主要 集中 在 北京 市 东南 部 、 天 


津 市 南部 及 河北 省 的 中 南部 , 低 低 聚 集 区 集中 分 布 在 西北 部 山区 。 研 究 区 高 低 聚 集 区 面积 均 呈 减 


小 趋势 ,不 显著 区 呈 扩 大 态势 。(3) 地 理 探测 器 结果 表明 不 同 区 域 的 主导 影 


响 因 子 不 同 ,北京 市 首 


A 
要 影响 因子 为 WOD 卫 ,其 次 为 人 口 密度 , 且 二 者 交互 作用 明显 。 天 津 市 主导 因子 为 风速 ,人 口 密度 、 
第 二 产业 生产 总 值 等 人 为 因子 的 作用 力也 较 大 ,风速 与 其 交互 作用 较 强 。 河 北 省 主导 因子 为 人 口 
密度 ,GDP 与 第 二 产业 生产 总 值 等 的 作用 力 次 之 ,整体 交互 作用 偏 弱 。 通 过 地 理 探测 器 解析 京 津 
ZW K AOD 空间 分 异 的 影响 机 理 , 其 结果 对 我 国 大 气 污染 治理 具有 重要 意义 。 
关 键 i: AOD; 时 空 变化 ; 地 理 探 测 器 ; 影响 因子 ; 京 津 费 地 区 


大 气 气 溶胶 是 指 大 量 悬 浮 在 空气 中 直径 小 于 
100 pm 的 液体 或 固体 微粒 , 它 不 仅 使 能 见 度 下 降 ， 
也 会 对 人 体 健 康 造 成 损害 ,是 大 气 污染 的 重要 组 成 
部 分 。 气 溶胶 光学 厚度 简称 AOD (Aerosol Opti- 
cal Depth ) 是 衡量 气 溶胶 的 重要 因子 ,可 以 用 来 表征 
大 气 浑 溃 度 或 气 溶 胺 的 总 含量 ,有 效 反 映 区 域 空气 
污染 程度 '"。 一 般 来 讲 ,40D 值 越 小 ,空气 中 气 溶 
胶 含 量 越 低 ,表明 大 气 越 清 洁 ; 反 之 ,40D 值 越 大 ， 
空气 中 气 溶胶 合 量 越 高 ,表明 大 气 越 混浊 。 目 前 
AOD 探测 主要 有 地 基 遥 感 与 卫星 和 遥感 两 种 手段 ,地 
面 设备 能 够 提供 气 溶胶 浓度 的 持续 变化 特征 ,但 不 


在 允许 误差 范围 内 ,可 以 被 用 于 全 球 范 围 的 气 溶胶 
变化 研究 ' 。 张 西 雅 等 5 、 孙 晓 雷 等 ”也 证 明 其 在 
京 津 旨 ,新疆 阿克苏 地 区 具有 较 好 的 适用 性 ,可 以 为 
区 域 尺 度 气 溶胶 研究 提供 更 精确 的 数据 支持 。 
SEF BEBE MX AOD 时 空 分 布 特征 及 其 影响 
机 理 已 经 有 了 不 少 研究 。 张 西 雅 等 中 对 京 津 蓉 地 
区 气 溶胶 时 空 分 布 状况 进行 研究 ,揭示 了 夜间 灯光 
数据 与 40D 具有 高 度 一 致 性 。 景 悦 等 ' PERE 
地 区 AOD 及 其 影响 因子 的 相关 分 析 表 明 京 津 蔓 地 
区 AOD 与 第 二 产业 生产 总 值 呈 明 显 正 相关 ,与 风速 
呈现 明显 负 相 关 。 刘 状 等 "对 中 国 北方 地 区 长 时 


能 得 到 良好 的 空间 覆盖 信息 ,难以 反映 其 变化 趋势 
的 空间 分 布 ,而 卫星 遥感 恰好 弥补 了 这 项 不 足 中 。 
近年 来 ,MODIS AOD 产品 在 大 气 污染 研究 领域 应 用 
最 为 广泛 ,其 标准 level 2 级 40D 产品 空间 分 辩 率 为 
10 km ,已 被 国内 外 大 量 学 者 证 明 具 有 显著 的 应 用 价 
值 ”-“ 2014 年 末 美 国 国家 航空 航天 局 发 布 了 更 
高 分 辩 率 的 3 km AOD 产品 ,该 产品 2/3 的 数据 都 
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期 的 气 溶胶 时 空 变化 特征 进行 分 析 , 表 明 人 口 密 度 
及 地 形 分 布 是 影响 AOD 空间 分 布 的 重要 原因 。 关 
于 气 溶胶 影响 因子 的 研究 ,通常 采用 相关 分 析 法 多 
元 线性 回归 等 方法 得 出 AOD 与 其 影响 因子 之 间 的 
线性 关系 “”" 。 而 地 理 探测 器 工具 对 变量 无 线性 
假设 ,可 以 探测 变量 与 影响 因子 之 间 的 所 有 关系 , 它 
是 基于 空间 分 异 理论 ,根据 要 素 的 空间 异 质 性 探测 


然 科 学 基金 青年 项 目 (41907194 ) ; R 


Lud 
=j 


和 自然 科学 基金 


作者 简介 : 景 悦 (1994 - +) , 女 , 硕 士 ,主要 研究 方向 为 大 气 环境 遥感 E-mail : 15222750696 @ 163. com 
通讯 作者 : 孙 艳 玲 (1977 - ) , 女 , 博 士 ,副教授 . E-mail; flyling99@ 163. com 


tó wy 


因子 解释 力 ,还 可 进一步 寻找 双 变 量 的 解释 能 力 , 其 
在 解释 变量 驱动 机 理 方面 具有 不 可 比拟 的 优势 。 目 
前 已 被 广泛 应 用 于 自然 .经 济 各 个 领域 55 -8 ,其 中 ， 
也 有 学 者 将 其 运用 到 空气 质量 4Q7 指数 的 风险 因子 
分 析 及 PM; ;驱动 因子 分 析 "” o (EE AOD 影响 
因子 识别 方面 研究 报道 较 少 ,仅见 郭 春 颖 利用 地 理 
探测 器 对 长 江 三 角 洲 地 区 AOD 的 风险 因子 进行 分 


析 ”" 。 因 此 ,本 文 利用 高 精度 气 溶胶 产品 MODO4_ 


3 上 数据, 结合 空间 自 相 关 模 型 与 地 理 探测 融 工 具 
科学 探测 和 诊断 了 京 津 副 地 区 气 溶 胶 时 空 规律 和 动 
力 机 制 ,有 效 的 验证 了 地 理 探 测 器 在 京 津 瘟 地 区 气 
浴 胶 驱动 机 理 方 面 的 适应 性 ,也 对 区 域 大 气 污 染 特 
征 认 识 与 治理 具有 重要 科学 意义 。 


1 数据 来 源 与 研究 方法 


1.1 数据 来 源 及 预 处 理 

(1) AOD; Terra MOD04 .3 上 数据 产品 是 从 美国 
国家 航空 航天 局 的 戈 达 德 宇 宙 飞 行 中 心 (https: // 
ladsweb. nascom. nasa. gov/ ) 网 站 上 下 载 得 到 ,时 间 


序列 为 2010 年 1 月 ~2016 年 12 月 ,空间 分 辨 率 为 
3 km, 时 间 分 辩 率 为 1 d, 波 段 为 550 nm 的 气 溶胶 


level 2 级 产品 ,除去 11 d 的 缺失 数据 , 共 得 到 2 547 
d 数 据 。40D 数据 利用 ENVI 的 IDL 编程 调用 
MCTK 工具 对 其 进行 投影 转换 批 处 理 , 再 进行 图 像 
拼接 HOMER ,在 剔除 无 效 值 及 0 值 之 后 生成 月 平 
均 数据 ,继而 得 到 年 平均 数据 。 

(2) 影响 因子 :分 为 自然 因子 和 人 为 因子 两 种 
类 型 。 其 中 ,自然 因子 包括 降水 量 、 相 对 湿度 .风速 、 
归 一 化 植被 指数 (NDVI) 。 气 象 数据 来 自 中 国 气 象 
数据 共享 网 (http : // data. cma. cn/ ) ,选择 克 里 金 插 
值 方法 对 风速 .降水 量 .相对 湿度 数据 进行 处 理 ”|。 
NDVI 数据 选用 MOD13A3 产品 数据 ,其 空间 分 辩 率 
为 1 km, 时 间 分 辩 率 为 15 d, 先 对 其 进行 投影 转换 、 
拼接 掩 膜 提取 等 预 处 理 ,进一步 生成 年 平均 数据 ， 
最 后 重 采样 生成 3 km 分 辨 率 的 栅 格 数据 。 人 为 因 
子 包括 国民 生产 总 值 (GDP) .第 二 产业 生产 总 值 和 
人 口 密度 。 数 据 来 自 于 北京 市 、 天 津 市 、 河 北 省 
2010—2017 年 统计 年 鉴 ,部 分 缺失 数据 从 各 省 统计 
局 官网 获得 ,统计 得 到 各 省 及 各 市 县 的 社会 经 济 数 
据 并 将 其 转化 成 矢量 面 数据 。 

PEAX AOD 值 及 各 影响 因子 的 7 a 平均 值 
空间 分 布 如 图 1 所 示 ,通过 地 图 判别 可 以 得 出 AOD 


及 各 影响 因子 空间 分 布 的 相似 性 与 差异 性 ,为 后 续 
进行 地 理 探测 的 结果 给 予 一 定 的 验证 。 
1.2 研究 方法 
1.2.1 室 间 自 相 关 分 析 地 理学 第 一 定律 表明 任 
何事 物 都 相关 ,只 是 相近 的 事物 关联 更 紧密 ,这 为 空 
间 自 相关 分 析 提 供 了 理论 基础 * i。 空间 自 相关 通 
篆 分 为 全 局 空间 自 相 关 和 局 部 空间 自 相 关 。 反 映 全 
局 空间 自 相 关 的 重要 指标 是 莫 兰 指数 (Moran 's I), 
Moran’s 1 >0 表示 空间 正 相 关 性 ,其 值 越 大 ,空间 相 
关 性 越 明 显 ,Moran's 7<0 表示 空间 负 相 关 性 ,其 值 
越 小 ,空间 差异 越 大 ,否则 ,Moran's7=0, 空 间 呈 随 
机 性 。 全 局 空间 自 相 关公 式 如 下 所 示 : 
W. ix = x)(x -x 
p Aane ee 


[Ly (x-7) TEE 


i=l j=l 


式 中 :n 为 单元 的 数量 ,为 AOD 栅 格 数据 转化 成 3 
km x3 km 的 矢量 面 数据 对 应 的 数量 ;xi 、w 为 单元 
ij AY AOD 的 年 均值 ; x 为 所 有 单元 的 40D 平均 值 ; 
W; 为 单元 i 的 空间 权重 矩阵 ,在 ArcGIS 中 用 原始 
方式 生成 默认 的 空间 权重 和 矩 了 泗 。 男 外 ,显著 性 水 平 
通过 Z 得 分 来 评定 ,其 计算 公式 为 : 
[1-E(D] 
VVar(7) 

=] (2) 
(n-1) 
Var(I) =E(P) - E(D)? 


AP: Z(1) ARKE Moran’s /的 显著 性 水 平 ;E(7) 
SN Moran’s 了 的 数学 期 望 ;Var(7) JH Moran’s 了 的 方 
差 。 

局 部 自 相 关 用 于 检测 要 素 的 高 值 或 低 值 的 局 部 
空间 聚集 性 ,一般 分 为 4 种 局 部 空间 相关 关系 :高 高 
聚集 表示 该 单元 与 其 周围 单元 的 AOD 值 均 相 对 较 
高 , 低 低 聚 集 指 该 单元 与 周围 区 域 的 40D 值 相对 较 
低 ,高 低 聚 集 指 高 值 中 心 被 低 值 所 包围 , 低 高 聚集 指 
低 值 中 心 被 高 值 所 包围 。 其 公式 如 下 所 示 : 


> yw, (x -¥) (x-7) 


peme (3) 
Ld (= -*) | 
式 中 :7 为 局 部 自 相 关 指 数 。 


Z(1) = 


E(I) = 
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1.2.2 地 理 探测 器 ”地理 探测 器 是 探测 要 素 的 空 
间 分 异性 并 揭示 其 背后 驱动 力 的 统计 工具 , 它 是 基 
于 空间 分 异 理论 ,将 自 变量 与 不 同 离散 分 类 的 自 变 
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一 定 


相似 性 "” 。 通 过 计算 q 统计 量 来 度量 要 素 的 空间 
异 质 性 以 及 探测 因子 解释 力 ,还 可 分 析 两 变量 真正 
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的 交互 作用 ,进一步 寻找 双 变量 的 解释 能 力 。9 值 


q% 的 了, 值 越 大 解释 能 
因子 探测 的 评价 指标 是 q 统计 量 , 其 公式 如 下 : 


n 
2 
È Nio; 
_ =l 


=p 
q No” 


值 域 范围 在 0 ~ 1 之 间 ,表示 自 变 量 工 解释 了 100 x 


就 越 强 。 


(4) 
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式 中 :=1,2,… 和 家 是 指 因子 蕊 的 分 类 ;N, 为 层 对 
应 的 单元 数 ;N 为 研究 区 总 单元 数目 ;o? 和 o? 分 别 
为 层 与 全 区 了 值 的 方差 。 了 代表 AOD ,驱动 因子 外 
包括 降雨 量 `. 相 对 湿度 .风速 . NDJ GDP, A O% 
E .第 二 产业 生产 总 值 ,用 3 km x3 km 空间 分 辩 率 
的 中 心 点 共计 24 223 个 网 格 点 将 变量 了 与 变量 
匹配 起 来 进行 探测 。 为 了 使 数据 在 空间 分 布 上 满足 
分 类 区 域 具有 显著 差异 性 ,根据 自然 断 点 法 .几何 间 
隔 法 以 及 百 分 位 数 法 对 影响 因子 进行 分 类 处 理 , 利 
用 地 理 探测 器 计算 每 一 种 分 类 的 g 值 ,最 终 选 择 解 
释 力 最 大 的 分 类 方法 及 组 数 作为 该 因子 最 佳 分 类 方 
A”, 

交互 探测 是 评估 X, 与 X 因子 相互 作用 是 否 会 
对 了 解释 力 产 生 影响 ,通过 分 别 比较 各 单 因 子 的 gq 
值 与 两 因子 有 加 后 的 g 值 ,可 以 判断 两 因子 交互 作 
用 的 强 弱 。 


2 结果 与 分 析 


2.1 京 津 费 地 区 AOD 的 时 间 变 化 特征 

通过 对 2010 一 2016 年 的 AOD 年 平均 数据 进行 
时 间 序 列 分 析 得 到 京 津 囊 地 区 ,北京 市 .天 津 市 、 河 
北 省 近 7 a AOD 的 变化 趋势 (图 2)。 结 果 表 明 
2010—2016 年 研究 区 整体 的 年 均 AOD 值 围绕 0. 83 
波动 ,7 a 间 AOD 值 呈 现 为 先 下 降 后 上 升 再 趋 于 稳 
定 的 状态 ,2010 年 40D 值 最 高 为 0.88,2012 年 AOD 
值 下 降 到 0.79 ,之 后 稳定 为 0.83。 北 京 市 的 7 a 
均 AOD 值 为 0.79 ,各 年 份 40D 值 均 低 于 研究 区 的 
整体 水 平 ,变化 趋势 与 研究 区 总 体 一 致 。 天 津 市 7 a 
平均 AOD 值 为 0.90 ,整体 40D 值 均 高 于 研究 区 的 
平均 水 平 ,呈现 “W” 型 变化 趋势 ,2010 年 AOD 值 最 
高 为 0.94 ,2012 年 出 现 第 一 个 谷 值 0. 87 ,2013 年 达 
到 峰值 0.90 ,随后 2014 年 达到 第 二 个 谷 值 0.88 ,之 
后 2a 又 呈 上 升 趋势 。 河 北 省 7 a 平均 AOD 值 为 
0.80 ,各 年 的 40D 值 均 小 于 研究 区 整体 水 平 ,40D 
值 在 2010 年 达到 最 高 为 0.84 ,2012 年 下 降 到 0. 76, 
之 后 在 2013 年 有 个 小 峰值 0. 80 ,随后 呈现 逐年 组 
慢 下 降 的 趋势 。 综 合 上 述 结果 可 以 看 出 ,各 地 区 
AOD 的 年 变化 趋势 与 研究 区 整体 大 致 相同 , 均 为 先 
下 降 后 上 升 再 小 幅 波 动 的 态势 。 天 津 市 的 40D 值 
为 研究 区 内 最 高 ,北京 市 最 低 , 河 北 省 的 40D 值 属 
于 中 间 水 平 。 
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图 2 2010—2016 年 京 津 费 地 区 AOD 的 时 间 变 化 特征 


Fig.2 Characteristics of the temporal variation of AOD 


in Beijing-Tianjin-Hebei region in 2010—2016 


2.2 HKEE AOD 的 空间 变化 特征 

全 局 自 相 关 分 析 结 果 表 明 ( 表 1) ,研究 时 段 内 
京 津 蔓 地 区 AOD 的 全 局 Moran's 了 指数 均 稳 定居 
高 ,其 中 北京 市 最 高 , 均 在 0.98 以 上 ,其 次 京 津 旨 及 
河北 省 地 区 为 0.96( +0.03), 天 津 市 相对 较 低 为 
0.95( +0.01) ,指示 各 地 区 AOD 的 空间 分 布 呈 现 显 
著 的 正 相关 ,同时 也 表明 京 津 束 地 区 AOD 呈现 集中 
连 片 的 状况 且 区 域 之 间 传 输 较 为 明显 
通过 对 2010—2016 FERE H X AOD 做 局 部 
自 相关 分 析 , 可 明确 研究 区 AOD 的 局 部 聚集 与 分 散 
情况 ,科学 解释 AOD 的 空间 变化 特征 。 根 据 结果 将 
显著 局 部 自 相关 的 区 域 划分 成 3 种 类 型 ,高 高 聚集 
型 . 低 低 聚集 型 及 不 显著 区 域 。 结 合 图 3 与 图 4 可 
知 , 各 地 区 AOD 的 局 部 聚集 特征 及 变化 趋势 不 尽 相 
同 。 北 京 市 的 低 低 聚集 型 面积 占 比 最 大 ,平均 占 北 
京 市 总 面积 的 45. 56% , 主要 分 布 在 西北 部 的 燕山 
山区 且 有 进一步 扩大 的 趋势 。 高 高 聚集 型 的 面积 次 
之 , 占 到 29.14% ,集中 分 布 在 城 6 区 及 大 兴 区 和 通 


o 


R1 2010—2016 年 京 津 费 地 区 AOD 全 局 自 相关 分 析 
Tab.1 Global spatial autocorrelation analysis of AOD in 


Beijing-Tianjin-Hebei region in 2010—2016 
全 局 自 相 关 ( Moran’s 7) 


年 份 z TER 
FEW ”北京 市 天 津 市 河北 省 

2010 0.93 0.98 0.95 0.93 
2011 0.97 0.98 0.94 0.97 
2012 0.95 0.99 0.94 0.95 
2013 0.96 0.99 0.96 0.96 
2014 0.95 0.99 0.95 0.96 
2015 0.96 0.99 0.96 0.95 
2016 0.98 0.99 0.94 0.98 
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最 小 为 25.29% ,主要 位 于 高 低 值 过 度 地 带 , 呈 减 小 
的 趋势 。 天 津 市 的 低 低 聚 集 型 平均 占 到 1.31% ,分 
布 区 域 极 少 但 呈现 增加 的 态势 。 高 高 聚集 型 平均 占 


图 3 2010 一 2016 年 京 津 豆 地 区 4OD 局 部 
自 相 关 分 析 


Fig.3 Local spatial autocorrelation analysis of AOD 
in Beijing—Tianjin—Hebei region in 2010—2016 


37.45% 的 不 显著 区 域 ,其 中 ,不 显著 区 域 有 进一步 
增加 的 可 能 ,而 高 高 聚集 型 呈 减 小 的 态势 。 总 体 而 
言 , 京 津 费 地 区 的 高 值 聚 集 区 、 低 值 聚集 区 及 不 显著 
区 域 各 占 总 面积 的 大 约 17/3 ,空间 上 大 致 呈现 东南 


天 津 市 总 面积 的 64. 01% ,分 布 在 武 清 一 宁 河 区 的 
以 南 地 区 , 呈 下 降 趋 势 。 不 显著 区 域 占 比 为 
34.67% , 主要 集中 在 北部 山区 且 无 明显 变化 趋势 。 
河北 省 的 低 低 聚 集 型 平均 占 比 为 33. 14% ,分 布 在 
太行 一 燕山 山地 地 区 ,有 逐渐 减 小 的 趋势 。 中 部 及 
东南 部 的 平原 地 区 分 布 有 27.51% 的 高 高 聚集 型 及 


高 西北 低 状况 。 北 京 市 东南 部 、 天 津 市 南部 .河北 省 
中 南部 始终 是 三 大 高 高 聚集 区 , 且 短 时 间 内 并 不 会 
发 生 改变 ,不 过 河北 省 南部 的 高 值 片 状 区 向 城市 中 
心 回 拢 .北京 市 天津 市 高 值 片 区 的 轻微 减少 可 能 

致 高 值 聚集 区 有 下 降 的 趋势 。 河 北 省 及 北京 市 的 西 
北部 低 低 聚集 区 的 减少 可 能 造成 京 津 豆 地 区 低 低 聚 
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图 4 2010—2016 年 京 津 蔓 地 区 局 部 聚集 区 面积 变化 趋势 


Fig.4 Area change trend of local gathering area in Beijing-Tianjin-Hebei region in 2010—2016 


集 区 面积 的 下 降 。 另 外 ,河北 省 南部 的 部 分 高 值 聚 
集 区 将 转化 为 不 显著 区 会 进一步 增加 京 津 豆 地 区 的 
不 显著 区 域 面积 。 
2.3 REKK AOD 的 影响 因子 分 析 

AOD 的 空间 变化 受气 修 、 地 形 、 人 类 活动 等 多 
种 因素 的 综合 影响 。 一 般 来 讲 , 相 对 湿度 增加 会 
使 气 溶胶 粒子 吸湿 增长 ,40D 增 大 ; 静 稳 天 气 的 环 
境 促进 污染 物 聚 集 , 大 风 天 气 会 加 快 扩散 ;降水 日 
数 较 多 利于 颗粒 物 的 吸收 与 沉降 ,40D 随 之 减 小 ; 
NDVI 反 映 地 表 植 被 覆盖 情况 ,植被 覆盖 好 对 调 
控 、 吸 收 大 气 污染 物 具有 明显 的 效果 ;GDP 反映 一 
个 地 区 经 济 发 展 程度 ,第 二 产业 生产 总 值 反映 工 
业 聚 集 情 况 ,这 些 都 与 40D 呈现 明显 的 正 相 关 ; 人 
口 密度 反映 地 区 人 口 活动 程度 ,人 口 越 多 ,表明 该 
地 区 人 们 生产 生活 频率 越 快 ,生活 的 废气 及 汽车 
尾气 排放 量 就 武大 。 因 此 ,研究 从 自然 和 人 为 两 
方面 选取 7 个 具有 代表 性 的 影响 因子 ,探测 其 对 
ERE HX. AOD 的 影响 机 理 。 利 用 地 理 探测 需 工 
有 具 首 先 对 不 同 地 区 各 影响 因子 贡献 率 进 行 分 析 
CPS) ,再 通过 因子 探测 与 交互 探测 得 出 京 津 豆 
地 区 不 同年 份 的 主导 驱动 因子 及 不 同 区 域 的 主导 
驱动 因子 ( 表 2) 。 


ri 


2.3.1 影响 因子 贡献 率 分 析 MARE K Ay 
响 因 子 的 贡献 率 可 以 看 出 (图 35) ,2010 一 2016 年 各 
因子 对 40D 平均 作用 强 弱 的 排名 依次 为 :人 口 密度 
(26.41% ) > 第 二 产业 生产 总 值 (22. 38% ) > GDP 
(21.82% ) > 相对 湿度 (12.71% ) > NDVI(7. 66% ) 
> 降水 量 (5.71% ) > 风速 (3.31% ) ,3 个 人 为 因子 
的 贡献 率 都 较 大 ,总 和 占 到 70. 61% ,说 明 人 为 因子 
是 该 区 域 40D 的 主要 驱动 力 。 从 时 间 序 列 来 看 ,各 
影响 因子 在 不 同年 份 对 40D 的 贡献 率 也 略 有 不 同 ， 
其 中 ,GDP 与 人 口 密度 在 各 年 的 贡献 率 基 本 变化 不 
大 。 第 二 产业 生产 总 值 在 2011 年 及 2014 年 贡献 率 
占 比 略微 下 降 到 18. 85% , 其余 年份 近似 为 
23.78% 。 降 水 量 对 AOD 的 贡献 率 在 2011 年 .2012 
年 较 高 为 11. 48% , 其 余年 份 都 影响 力 较 低 , 为 
3.40% 。 相 对 湿度 贡献 率 的 变化 呈现 倒 “U” 型 趋 
势 ,总体 平 均 贡 献 率 为 12.71% ,最 高 贡献 率 出 现在 
2013 4654 18.72% 。NDVI 的 变化 情况 与 降水 量 呈 
相反 态势 ,NDVI 的 贡献 率 恰好 在 2011 和 2012 年 降 
到 最 低 4.50% ,其 余 5 a 平 均 贡 献 率 为 8.91% 。 风 
EXT AOD 贡献 率 的 变化 趋势 与 降水 量 相似 ,整体 呈 
现 单 峰 分 布 ,峰值 出 现在 2010 年 ,贡献 率 为 6.22% 。 

北京 市 驱动 因子 的 探测 结果 显示 (图 5) ,各 因 
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子平 均 贡 献 率 强 弱 依次 为 :VDIT(21.32% ) > AO 
密度 (17.72% ) > GDP(15.76% ) > 第 二 产业 生产 
总 值 (14. 39% ) > 风速 (13. 43% ) > 相对 湿度 
(10.05% ) > 降水 量 (7.33% ) ,前 三 者 贡献 率 总 和 
占 比 为 54. 80% ,说 明 北 京 市 的 NDVI 与 其 AOD 变 
化 规律 较 一 致 ,日 人 口 高 密度 区 及 经 济 快速 发 展区 
的 地 理 分 布 与 40D 高 值 区 具有 高 度 的 空间 耦合 
从 不 同年 份 各 因子 贡献 率 的 变化 中 看 出 , NDVI A 
口 密度 GDP 相对 比较 稳定 。 第 二 产业 生产 总 值 的 
贡献 率 7 a 间 出 现 小 幅 震 荡 , 其 中 2012 年 .2014 年 、 
2015 年 出 现 了 较 低 贡献 率 为 12. 65% ( +0.76), 其 
余年 份 贡献 率 的 均值 达到 15. 69% 。 降 水 量 的 贡献 
率 在 2010 一 2015 年 均 稳定 在 8.33% ,到 2016 年 急 
剧 下 降 到 1. 25% 。 相 对 湿度 的 贡献 率 在 2010 一 
2011 年 较 低 ,为 5.01% ,在 2012 年 达到 一 个 峰值 ， 
贡献 率 高 达 16. 88% ,但 随后 又 下 降 并 稳定 在 
10. 87% ( +2.54)。 风 速 对 AOD 的 贡献 率 在 2013 
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年 有 一 个 谷 值 为 9. 31% , 其 余年 份 的 贡献 率 约 为 
14.72% 。 

天 津 市 因子 探测 的 结果 具有 自身 的 独特 性 (图 
5),7 a 平均 的 各 因子 贡献 率 排 名 依次 为 风速 
(21.24%) > 人 口 密度 (19. 05%) > NDVI 
(16.73% ) >GDP(13. 56% ) > 第 二 产业 生产 总 值 
(12. 88% ) > 降水 量 (9. 16% ) > 相对 湿度 (7.38% ) 。 
说 明天 津 市 人 口 高 密度 区 及 强风 速 地 区 的 40D 值 
较 高 ,而 北部 人 口 稀 玖 、 山 区 植被 覆盖 较 好 恰好 与 
AOD 低 值 区 的 范围 近似 。 从 时 间 序 列 来 看 ,风速 、 
NDVI 人 人口 密度 3 个 因子 的 贡献 率 变 化 不 大 ,较为 
稳定 。GDP 的 贡献 率 在 7 a 间 有 小 幅 震 荡 ,2010 年 
及 2012 年 贡献 率 偏 低 为 10.49% ,其 余年 份 平均 贡 
献 率 基 本 一 致 为 14.79% 。 第 二 产业 生产 总 值 的 贡 
献 率 呈 U 型 分 布 ,在 2013 年 .2014 年 有 一 个 低谷 值 
为 8.66% ,其 他 年 份 平 均 为 14.57% 。 降 水 量 的 变 
化 趋势 波动 频繁 ,2011 年 与 2015 年 降水 量 的 贡献 
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Fig.5 Contribution rate of influence factors of AOD in Beijing-Tianjin-Hebei region in 2010—2016 
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率 仅 为 1.57( +0.48) ,几乎 没有 影响 ,2016 年 贡献 
力 轻微 上 升 到 3.99% ,其 余年 份 贡献 率 较 高 稳定 在 
14.25% 。 相 对 湿度 对 天 津 市 40D 的 影响 总 体 呈 现 
逐渐 上 升 的 过 程 ,2016 年 贡献 率 达 到 最 高 为 
11.88% 。 

河北 省 单 因子 贡献 率 的 分 析 中 (图 5) ,各 因子 
7 a 平 均 的 影响 力 排 序 依次 为 人 口 密 度 (25.38% ) 
> GDP(23. 47% ) > 第 二 产业 生产 总 值 (21.75% ) 
> 相对 湿度 (13. 57% ) > 降水 量 (6.56% ) > NDVI 
(4.77% ) > 风速 (4. 50% ) 。 河 北 省 的 因子 探测 结 
果 也 显示 人 为 因素 三 者 总 体 贡 献 率 达 到 70. 06% ， 
说 明 经 济 发 展 活动 及 人 类 生活 排放 是 影响 河北 省 
AOD 的 主要 决定 因素 。 不 同年 份 各 影响 因子 的 变 
化 也 各 不 相同 ,人 口 密 度 与 GDP 两 个 人 为 因子 的 逐 
年 变化 最 为 稳定 ,而 第 二 产业 生产 总 值 的 贡献 率 在 
2012 年 出 现 一 个 低谷 值 为 15. 61% ,其 余年 份 均 较 
高 为 22.77% 。 降 水 量 贡 献 率 的 变化 趋势 呈 先 上 升 
后 下 降 的 过 程 ,2012 年 达到 峰值 为 11. 99% ,2016 
年 降 到 最 低 值 为 2.78% 。 相 对 湿度 的 时 波动 上 升 
的 趋势 , 2016 年 达到 贡献 率 的 最 高 值 17. 81% 。 
NDVI 的 变化 趋势 呈现 轻微 震荡 ,2010 年 与 2015 年 
NDVI 的 贡献 率 较 高 为 7.06( +0.35) ,其余 年份 的 
贡献 率 都 较 低 ,平均 为 3.86% 。 风 速 的 变化 趋势 与 
降水 量 一 致 ,不 过 ,其 整体 平均 贡献 率 都 较 低 ,最 高 
才 达 到 7.43%。 

河北 省 因子 探测 结果 与 京 津 蔓 整 体 结果 类 似 ， 
影响 因子 的 逐年 变化 规律 也 大 致 相同 。 而 北京 市 、 
天 津 市 各 因子 贡献 率 确 大 有 不 同 ,证 明了 全 局 因子 
探测 与 局 部 因子 探测 的 重要 性 ,同时 也 启示 区 域 大 
气 污染 治理 应 该 因地制宜 . 找 准 各 区 问题 的 关键 。 
2.3.2 主导 因子 分 析 MEAW AOD 的 成 因 复 
杂 、 驱 动因 子 多 样 , 准 确 把 握 主 导 因 子 是 解决 区 域 气 
溶胶 污染 的 关键 。 将 地 理 探测 器 分 别 运用 于 京 津 蔓 
地 区 ,北京 市 .天津 市 河北 省 4 个 区 域 ,探测 各 个 地 
区 逐年 的 主导 因子 及 交互 主导 因子 ( 表 2) ,结果 表 
明 ,最 大 9 值 对 应 的 因子 随 着 不 同年 份 不同 地 区 会 
产生 一 定 的 变化 ,其 中 不 同 地 区 之 间 的 主导 因子 差 
异 巨 大 。 就 单一 主导 因子 而 言 , 京 津 间 地 区 2010 
年 .2011 年 .2014 一 2016 年 主导 因子 为 人 口 密度 ， 
2012 年 和 2013 年 主导 因子 为 第 二 产业 生产 总 值 ,7 
a 间 单一 主导 因子 的 平均 q 值 为 0.39, 其 中 2010 
年 .2011 年 .2016 年 最 大 4 值 接近 0.5 ,说 明 该 时 期 


表 2 2010—2016 年 京 津 费 地 区 主导 因子 及 交互 主导 
因子 影响 强度 (4 值 ) 表 
Tab.2 Influence of dominant factors and interacting 
dominant factors in Beijing-Tianjin-Hebei region in 


2010—2016(q value) 


m 2010 ”人口 密度 0.45 人 口 密度 mmNDT7 0. 62 
地 区 2011 人口 密度 0.48 人口 密 度 门 风速 0. 56 
2012 0.37 oe 0.46 
第 二 产 第 二 产业 生产 
2013 or 0.32 ae ESEN 0.43 
2014 ”人 口 密度 0.31 AFAR BEM KE 0. 44 
2015 ”人 口 密度 0.32 AODREN KE 0. 42 
2016 AHR 0.50 ”人 口 密度 六 相对 湿度 0.63 
北京 市 2010 NDVI 0.63 NDVIN 风速 0.78 
2011 NDVI 0.65 ”NDVIN 人 口 密度 0. 84 
2012 NDVI 0.64 NDVIN 人 口 密度 0. 83 
2013 NDVI 0.67 NDVIN 人 口 密度 0. 82 
2014 NDVI 0.67 NDVIN 人 口 密度 0. 85 
2015 NDVI 0.64 NDVIN 人 口 密 度 0. 84 
2016 NDVI 0.72 NDVIN 人 口 密度 0. 85 
天 津 市 2010 ”风速 0.40 ”风速 人口 密 度 0. 46 
2011 ”风速 0.39 ”风速 但 相对 湿度 0.50 
2012 ”风速 0.46 ”风速 人 人 口 密 度 0. 54 
2013 ”风速 0. 45 速 儿 相对 湿度 0. 68 
2014 ”风速 0.35 ”风速 儿 相 对 湿度 0. 52 
2015 ”风速 0.36 ee 0.51 
2016 “人口 密 度 0.51 ”风速 站 人 口 密度 0. 69 
河北 省 2010 — ai oe eee 
2011 人口 密度 0.48 ”人 口 密度 Nn 降水 量 0.54 
2012 ”人口 密度 0.35 A aE 风速 0. 43 
2013 “人口 密度 0.30 人口 密度 站 相对 湿度 0. 40 
2014 ” 人口 密 度 0.30 ”人 口 密度 Nn 降 水 量 0.44 
2015 A FABRE 0.32 AHE NKE 0. 40 
2016 “人口 密 度 0.52 人口 密度 相对 湿度 0.64 
对 应 的 单一 主导 因子 可 以 解释 了 50% 的 40D 态势 。 


北京 市 的 单一 主导 因子 7 a 均 为 NDVI, PH q 值 高 
达 0.66, 各 年 份 g 值 均 大 于 0.6, 最 高 9 值 为 2016 
年 的 0.72 ,表明 2016 年 的 NDVI 分 布 状 况 可 以 反映 
72% 的 40D 空间 规律 。 天 津 市 2010 一 2015 年 均 为 
风速 作用 强度 最 大 , BI q 值 变化 不 大 ,稳定 在 
0.40( +0.06) 上 下 ,2016 年 主导 因子 转变 为 人 口 密 
度 ,gq 值 达 到 0.51 ,可 以 看 出 该 年 份 人 为 活动 影响 已 
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经 超过 了 风速 对 40D 的 作用 强度 。 河 北 省 2010 年 
的 主导 因子 为 第 二 产业 生产 总 值 ,2011 年 之 后 均 为 
人 口 密度 ,7 a 平均 g 值 在 0.35 上 下 浮动 , 较 大 7 值 
出 现在 2011 年 .2016 年 ,这 2 a 主导 因子 可 以 解释 
接近 一 半 的 AOD 分 布 情况 。 在 交互 主导 因子 中 ,最 
K q 值 的 交互 主导 因子 都 是 该 年 度 的 单一 主导 因子 
与 另 一 因子 交互 的 结果 ,而 且 任 意 两 个 因子 交互 后 
的 作用 强度 都 明显 大 于 单一 因子 ,但 不 同 地 区 交互 
作用 强 弱 略 有 差异 , 京 津 蔓 地 区 主导 交互 因子 的 平 
I q 值 为 0.51 , 比 单一 主导 因子 增加 了 0.12 。 北 京 
市 主导 交互 因子 的 q 值 增加 幅度 较 大 ,交互 后 的 平 
I q 值 为 0. 83 。 天 津 市 主导 因子 交互 作用 后 9 值 增 
加 了 0. 14 ,平均 达到 56% 的 解释 能 力 。 河 北 省 交互 
作用 最 不 显著 ,平均 q 值 增加 了 0. 09 ,交互 后 平均 
解释 能 力 为 44% 。 可 以 看 出 ,交互 作用 在 北京 市 最 
为 显著 , 京 津 蔓 地 区 及 天 津 市 的 交互 作用 次 之 ,河北 
省 的 作用 最 不 明显 ,因此 ,今后 对 于 北京 市 大 气 污染 
的 防 控 更 应 该 注意 NDVI 与 人 口 密度 的 交互 作用 ， 
天 津 市 应 注重 风速 与 人 口 密度 .相对 湿度 .第 二 产业 
等 的 复合 作用 力 。 

综合 上 述 不 同 地 区 的 因子 探测 的 结果 ,进一步 
提出 各 个 区 域 精准 治理 大 气 污染 的 对 策 与 建议 。 其 


3 讨论 


本 文采 用 AOD 遥感 数据 来 反映 区 域 空气 污染 
情况 而 未 采用 AOD 反 演 的 PM, ;数据 进行 研究 , 主 
要 原因 在 于 AOD 本 身 就 可 有 效 反映 区 域 气 溶胶 ( 直 
径 小 于 100 um 的 液体 或 固体 微粒 ) 状况 ,而且 精 度 
较 高 "| 。 通 过 AOD 反 演 PM, ;后 ,多 种 不 确定 性 因 
素 会 导致 反 演 的 PM, ;浓度 存在 一 定 的 误差 ,目前 比 
较 先 进 的 模型 ( 时空 地理 加 权 模 型 ,两 步 又 混合 效 
应 模型 . 贝 叶 斯 模型 等 ) 反 演 PM, ;浓度 的 平均 预测 
误差 也 达到 12 ~ 18 pg m ,其 至 不 同 模型 
添加 不 同 的 辅助 因子 ,会 得 到 不 同 的 结果 。 同 时 李 
成 才 等 ”也 证 明了 高 分 辩 率 的 MODIS AOD 数据 研 
究 城市 地 区 的 空气 污染 具有 显著 的 可 行 性 。 因 此 本 
研究 选择 NASA 提供 的 MODIS AOD 数据 直接 分 析 。 

另外 ,40D 定义 为 介质 的 消光 系数 在 垂直 方向 
上 的 积分 ,是 描述 一 定 的 垂直 高 度 上 和气 溶胶 对 光 的 
削减 作用 。 若 将 AOD 垂直 订正 为 地 面 消光 系数 ,与 
地 面 气象 因子 及 人 为 因子 进行 影响 分 析 显 然 更 为 合 
理 , 但 由 于 获取 2 547 d 不 同 地 区 的 边界 层 高 度数 据 
困难 较 大 ,本 文 利用 整个 区 域 的 年 平均 边界 层 高 度 
数据 进行 试验 ,研究 结果 显示 京 津 蔓 地 区 年 际 边界 


中 ,北京 市 与 天 津 市 属于 自然 一 人 为 因子 双重 驱动 
型 ,河北 省 属于 人 为 因子 驱动 型 。 北 京 市 首要 主导 
因子 为 NDVI, 其 次 为 人 口 密 度 , 该 因子 对 AOD 的 影 
响 力 稳定 居 高 且 与 NDVI 交互 作用 巨大 。 启 示 北 京 
市 应 该 提高 植被 覆盖 率 , 有 效 发 挥 其 对 气 涂 胶 粒子 
的 稀释 作用 ,降低 污染 ,并 控制 人 类 活动 的 影响 , 推 
崇 使 用 公共 交通 ,减少 汽车 尾气 的 排放 。 天 津 市 的 
主导 因子 为 风速 ,人 口 密度 及 NDVI 的 影响 力 次 之 ， 
GDP 的 作用 力 逐 年 升 高 , 且 风 速 与 人 为 因子 、 相 对 
湿度 的 复合 作用 力 显 车。 进一步 提出 天 津 市 应 有 效 
利用 风力 这 一 自然 条 件 ,控制 污染 企业 的 合理 布局 ， 
促进 本 地 污染 物 扩散 ,增加 植被 覆盖 度 有 效 防治 扬 
尘 地 块 ,减少 风蚀 起 尘 ,最 后 应 重点 控制 人 类 生产 生 
活 废气 的 排放 ,防止 人 为 污染 的 进一步 加 剧 。 河 北 
省 属于 人 为 因子 驱动 型 ,人 口 密度 是 AOD 的 首要 影 
响 因 子 ,GDP 及 第 二 产业 的 作用 力 次 之 ,相对 湿度 


层 高 度 变化 不 大 ,用 其 进行 垂直 订正 的 效果 不 明显 ， 
最 终 将 其 统一 视 为 AOD 估算 的 系统 误差 。 在 未 来 
的 研究 中 ,应 该 使 用 每 日 不 同 地 区 的 边界 层 高 度数 
据 进 行 垂直 订正 ,才能 保证 订正 结果 的 准确 性 。 

与 前 人 对 京 津 费 地 区 AOD 的 影响 因子 的 研究 
结果 对 比 可 知 ,本文 结 果 表 明 京 津 蔓 全 区 及 河北 省 
的 首要 影响 因子 为 人 口 密度 ,这 与 张 西 雅 等 |、 刘 
状 等 中 "研究 结果 一 致 。 但 本 文 基于 不 同 区 域 的 探 
测 结果 显示 北京 和 天 津 的 自然 因子 也 起 到 较 大 的 作 
用 ,这 与 刘 壮 等 认为 40D 受 自 然 源 影响 较 小 的 结论 
不 一 致 ,造成 两 者 结论 不 同 的 原因 可 能 是 刘 状 等 并 
未 进行 40D 与 自然 因子 的 定量 分 析 , 只 是 认为 人 口 
密度 增长 与 40D 增长 趋势 及 范围 十 分 一 致 ,进一步 
推测 说 明 自 然 因子 影响 较 小 ,而 本 文 对 自然 与 人 为 
因子 均 进 行 定量 分 析 ,结果 相对 更 为 可 靠 。 

AOD 与 其 影响 因子 的 Pearson 相关 分 析 结 果 显 


的 影响 力 有 逐渐 增加 的 趋势 ,针对 上 述 主导 因素 , 河 
北 省 应 着 重 降低 人 类 活动 对 污染 物 的 排放 ,杜绝 散 
煤 烟 煤 的 使 用 ,发 展 清洁 能 源 ,提倡 公众 绿色 生活 ， 
同时 应 该 调整 产业 结构 ,促进 第 二 产业 的 能 源 优化 ， 
鼓励 发 展 第 三 产业 等 。 


示 , 京 津 蔓 地 区 及 河北 省 是 人 为 因子 占 主导 ,与 地 理 
探测 结果 一 致 ,而 天 津 市 与 北京 市 相关 性 较 明显 的 
均 为 自然 因子 ,地 理 探测 器 的 结果 则 表明 两 地 区 
AOD 的 空间 分 异 是 受 自 然 与 人 为 双重 作用 的 结果 。 
两 者 结果 既 有 相似 性 也 存在 一 定 的 差异 性 ,这 主要 
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与 两 个 统计 方法 的 原理 不 同 有 关 ,Pearson 相关 系数 
分 析 法 主要 是 衡量 AOD 及 影响 因子 之 间 线 性 关系 ， 
而 地 理 探测 器 则 对 变量 无 线性 假设 ,可 以 探测 变量 
与 影响 因子 之 间 的 所 有 关系 。 且 京 津 费 地 区 AOD 
及 其 影响 因子 的 空间 差异 较 大 ,将 40D 与 其 影响 
子 看 成 一 个 整体 进行 Pearson 相关 分 析 , 显然 忽视 
了 要 素 的 空间 异 质 性 , 而 地 理 探 测 需 基 于 空间 分 异 
理论 ,探测 导致 要 素 空间 分 异 的 驱动 因子 ,从 这 个 角 
度 看 ,地 理 探测 器 对 AOD 的 影响 机 理 探 测 可 能 更 加 
合理 。 目 前 对 京 津 蔓 地 区 AOD 驱动 因子 的 空间 探 
测 前 所 未 有 ,因此 ,本文 的 研究 结果 对 该 地 区 的 大 气 
环境 监测 以 及 工农 业 生产 布局 等 具有 重要 意义 。 


4 结论 


(1) 2010—2016 EREHE X AOD 的 年 均值 
为 0. 83 ,其 中 ,天 津 市 40D 的 年 均值 为 研究 区 最 高 ， 
河北 省 次 之 ,北京 市 最 低 。 研 究 区 整体 及 各 分 区 
AOD 的 年 变化 趋势 大 体 一 致 ,呈现 先 下 降 后 上 升 再 
小 幅 波 动 的 态势 。 

(2) 空间 自 相 关 结 果 表 明 研 究 时 段 内 京 津 蔓 地 
区 AOD 空间 分 布 呈 现 显著 的 正 相 关 。 局 部 高 高 聚 
RX \ 低 低 聚 集 区 及 不 显著 区 域 各 占 总 面积 的 大 约 
1/3。 北 京 市 东南 部 、 天 津 市 南部 、 河 北 省 中 南部 始 
终 是 三 大 高 高 聚集 区 ,呈现 轻微 缩小 的 态势 , 低 低 聚 
集 区 位 于 西北 部 山区 ,呈现 减 小 趋势 ,其 余地 区 均 为 
不 显著 区 域 , 面 积 呈 扩 大 态势 。 

(3) 地 理 探测 器 的 结果 显示 北京 市 首要 影响 
子 为 NDVI ,其 次 为 人 口 密度 , 且 二 者 交互 作用 后 对 
AOD 影响 巨大 。 天 津 市 主导 因子 为 风速 ,人 口 密 
度 .第 二 产业 生产 总 值 等 人 为 因子 的 作用 力 逐 渐 增 
大 ,风速 与 其 交互 作用 较 强 。 河 北 省 主导 因子 为 人 
口 密度 , GDP 与 第 二 产业 生产 总 值 等 的 作用 力 次 
之 ,整体 交互 作用 偏 弱 。 启 示 北 京 市 与 天 津 市 应 有 
效 发 挥 自然 因子 的 正 效 应 ,并 同时 防 控 好 人 类 活动 
的 影响 ,河北 省 应 重点 减少 人 为 活动 的 排放 。 
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Abstract; The Beijing-Tianjin-Hebei region , which is located in the Bohai Rim, is the most economically devel- 
oped region in northern China. It is also one of the most polluted regions in China. Aerosol optical depth (AOD) can 
effectively reflect the degree of air pollution. In this study ,the MODIS 3 km AOD remote sensing data was analyzed 
and the research results would be of great significance to China’s air pollution. Firstly ,the temporal variation charac- 
teristics and spatial autocorrelation characteristics of AOD in Beijing-Tianjin-Hebei region from 2010 to 2016 were 
analyzed. The factors of precipitation, wind speed, relative humidity, normalized difference vegetation index ( ND- 
VI) ,gross domestic product (GDP) ,secondary industry GDP, population density were selected as the affecting fac- 
tors of AOD, using geographical detection to analyze the trend of the factors contribution ratio and its dominant af- 
fecting factors in different regions and different times in Beijing-Tianjin-Hebei region. The results showed as follows; 
(1) The annual average AOD of the Beijing-Tianjin-Hebei region was 0. 83 from 2010 to 2016. The annual average 
AOD of Tianjin was the highest in the study area, followed by Hebei and Beijing. The annual variation trend of AOD 
was shown an overall state of decline first followed by rise and then a small fluctuation change and a similar status 
in each local region as well. (2) The spatial autocorrelation analysis indicated a significant positive correlation a- 
mong spatial distribution of AOD in Beijing-Tianjin-Hebei region. Local high-accumulation areas were mainly con- 
centrated in the southeastern part of Beijing ,southern part of Tianjin, and central-southern part of Hebei Province. 
The low-accumulation areas were concentrated in the mountains of the northwest. The area of the high and low accu- 
mulation areas in Beijing-Tianjin-Hebei region showed a decreasing trend, while the area of the non-significant areas 
showed an increasing trend. (3) Geographic detection analysis indicated that the primary influencing factor was 
NDVI in Beijing, followed by population density ,and their interaction was significant. The dominant factor was wind 
speed in Tianjin. The effect of human factors such as population density and the GDP of the secondary industry were 
also important. The effect of interactions between wind speed and the above factors were also significant. The domi- 
nant factor was population density in Hebei Province, followed by GDP and the secondary industry’s GDP, and the 
overall interaction was relatively weak. In the past, the correlation analysis of the influencing factors was based on 
linear relationship of the overall area, while the geographical detection not only has wireless assumptions, but also 
reveals the linear, nonlinear and spatial relationships of the driving factors. From this perspective , the geographical 
detection is more reliable in studying the influencing factors. This study applied geographical detection to analysis of 
the impacting factors of AOD in Beijing-Tianjin-Hebei region , which proved its feasibility in this field and this re- 
gion. At the same time, the results also had important reference value for air pollution prevention, industrial and ag- 
ricultural layout and urban construction in Beijing-Tianjin-Hebei region. 
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